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ABSTRACT 

Pulsed 7-ray emission from millisecond pulsars (MSPs) has been detected by the sensitive Fermi, which 
sheds light on studies of the emission region and mechanism. In particular, the specific patterns of radio and 
7-ray emission from PSR J0101-6422 challenge the popular pulsar models, e.g. outer gap and two-pole caustic 
(TPC) models. Using the three dimension (3D) annular gap model, we have jointly simulated radio and 7- 
ray light curves for three representative MSPs (PSR J0034-0534, PSR J0101-6422 and PSR J0437-4715) with 
distinct radio phase lags and present the best simulated results for these MSPs, particularly for PSR J0101-6422 
with complex radio and 7-ray pulse profiles and for PSR J0437-4715 with a radio interpulse. It is found that 
both the 7-ray and radio emission originate from the annular gap region located in only one magnetic pole, and 
the radio emission region is not primarily lower than the 7-ray region in most cases. In addition, the annular 
gap model with a small magnetic inclination angle instead of "orthogonal rotator" can account for MSPs' radio 
interpulse with a large phase separation from the main pulse. The annular gap model is a self-consistent model 
not only for young pulsars but also MSPs, and multi-wavelength light curves can be fundamentally explained 
by this model. 

Subject headings: pulsars: general — gamma rays: stars — pulsars: individual (J0034-0534, J0101-6422, 
J0437-4715) — radiation mechanisms: non-thermal — acceleration of particles 



1. INTRODUCTION 

MSPs are a population of old neutron stars with short spin 
period Pq < 20 ms and small characteristic magnetic field 



Bo < 10 G (actually small period derivative P < 10 s s 



They are believed to be for med from the recycling (accretion 
spin-u p ) process in binaries dBhattacharva & van den Heuvell 
11991b . Furthermore, there is another possible formation 
channel for MSPs (even for sub-millisecond pulsars), that 
is accretion-induced collapse of white dwarfs in a binary 
(iNomoto & Kondoll99UlDu et al.l2009l Because of the high 
stability of their rotations, MSPs have great potential for ap- 
plication: autonomous deep-space navigation, pulsar-based 
time scale, low-frequency gravitational wave detection and so 
on. 

When the first MSP, PSR B 1 937+21, ha d been discov- 
ered in 1982 dBacker et al.lll982l) . IUsovl(ll983h soon predicted 
that this pulsar could emi t 7-rays via synchrotron radiation 
on the order of 100 GeV. iBhattacharva & Srinivasanl (11991b 
subsequently calculated the 7-ray luminosity of MSPs, esti- 
mated their contribution to the diffuse 7-ray background of 
the Milky Way, and finally discussed detectability of MSPs as 
point 7-ray s ources. Using the d ata of Energetic Gamma Ray 
Experiment, iKuiper et al.l (120001) showed circumstantial evi- 
dence for the likely detection of pulsed 7-ray emission from 
a MSP, PSR J02 18+4232, which was regarded as a 7-ray pul- 
sar candidate for a long time. MSPs' 7-ray emission was not 
observationally confirmed until Fermi with a sensitive Large 
Area Telescope (LAT) launched in June of 2008. Using eight- 
month data of Fermi LAT, eigh t 7-ray MSPs have already 
been detected dAbdo et al. 2009). This number has grown to 
more than 40 up to date dGuillemot et al. |[2012al) . and still 
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increases. 

From observations, MSPs are analogous to young pulsars, 
which have multi-wavelength pulsed emission from radio to 
7-ray band. Do MSPs and young pulsars share a simple model 
that contains similar emission region and acceleration mecha- 
nism to self-consistently explain their multi-wavelength emis- 
sion? More and more multi-wavelength data with high preci- 
sion give us opportunities to obtain remarkable insights of the 
magnetospheric physics. The multi-wavelength study is a key 
method to discriminate the various pulsar non-thermal emis- 
sion models for both MSPs and young pulsars. 

Initially aiming to explain the high-energy pulsed 
emission from young pulsars, four traditional magne- 
tospheric gap models have been suggested to study 
pulsed high energy emission o f puls ars: the polar 
cap mo del dbaughertv & Harding] 119941), the out e r gap 
model (ICheng et al.l Il9 86t iRoman i & Yadigaroglul 119951: 
Zhang & Chendl 19971). TPC /slot gap model dDvk s & Rudakl 
20031: lHarding et all |2008[). and the annular gap model 
(lOiao etldl 12004 120071: iDu et al.ll2010l) . To distinguish these 
pulsar models, the most important issues are the induce- 
ments of acceleration electric field region and relat ed emis- 
sion mecha nisms to emanate high-energy photons dDu et al.l 
|2011L|20T2|) . One of the key discrepancies of the mentioned 
emission models is one-pole or two-pole emission pattern 
which depends on two corresponding geometry parameters: 
m agnetic incl i nation angle a and viewing angle £. 

iBulik et all ( 120001) adopted the polar-cap model to calculate 
the 7-ray emission from MSPs. They pointed out that cur- 
vature radiation of primary particles contributed to the MeV- 
to-GeV band, while the synchrotr on radiation ar i sing f rom 
pairs dominated only below 1 Me V. lHarding et al.l (120051) de- 
veloped the pair-starved polar cap model and obtain similar 
spectral conclusion for high-energy emission from MSPs as 
above. In this model the accelerating field is not screened and 
the entire open volume is avai lable for particle accele ration 
and emission of gamma rays. iZhang & Chengl d200l used 



2 



Duetal. 2012 



the outer gap model with multi-pole magnetic field to model 
the X-ray and 7-ray spectra for four MSPs, an d the predicted 
result s basically agree with the observations dHarding et alj 
[2005b. 

Along with radio observations supplying us with excellent 
timing solutions for Fermi MSPs, the derived 7-ray and ra- 
dio light curves with high signal-to-noise allow us to do joint 
simulat ion which ca n justif y the pulsar emission models. Re- 
cently, IVenter et al.l (120091) simulated both radio and 7-ray 
light curves for MSPs in the pair-starved polar cap, TPC and 
outer gap models, and they found that most of their simulated 
light curves are well explained by the TPC and outer gap mod- 
els. They especially simulated light curves for a minor group 
of MSPs with pha se-aligned radi o and 7-ray p u lse pr ofiles 
( IVenter et al.l 120121) . In addition, Uohnson et at] d2012l) also 
used the geometric slot gap, outer gap model or pair-starved 
polar cap model to fit 7-ray and radio light curves for three 
MSPs. 

A MSP PSR J0101-6422 with complex radio and 7-ray 
light curves, chal lenge the popular TPC and outer gap models 
dKerr et al.ll2012l) . It is found that neither of the two models 
can faithfully reproduce the observed light curves and phase 
lags. For such a complex radio profile the simple beam model 
they used may have been insufficient. 

In this paper, we use a 3D annular gap model to study both 
radio and 7-ray light curves for three MSPs which stand for 
the relevant types of MSPs with different radio lags. In § 2, 
the annular gap and core gap will be simply introduced and 
the acceleration electric field in the annular gap is calculated. 
In § 3, we jointly simulate radio and 7-ray band light curves 
for PSR J0034-0534, J0101-6422 and J0437-4715. The radio 
phase lags are identified and the reasons for them explored. 
The conclusions and discussions are shown in § 4. 

2. THE ANNULAR GAP AND CORE GAP FOR MSPS 

2.1. The Definition of Annular Gap and Core Gap 
In a pulsar magnetospher e dGoldreich & Juliarj_ 119691: 



IRuderman & Sutherland|[T975]) . the critical field lined divide 
a polar cap into two distinct parts: the annular gap regi on and 
the core gap region (See Figure 1 of iDu et alj 120121) . The 
annular gap is constrained between the critical and last open 
field lines, and the core gap is around th e magnetic axis and 
within critical field lines dDu et al.ll20lol) . The size of the po- 
lar cap decreases with increasing spin period, but can be quite 
large for MSPs and young pulsars. The annular gap width 
is correspondingly larger for short spin-period pulsars, and 
it varies with the magnetic azimuth ip. For an anti-parallel 
rotator and ip = 0°, the radii of the annular polar region is 
farm = r pc -r cole = 0.26R(SIR / c) x l 2 dDu et al.H2010l) . here R is a 
pulsar's radius, and ft is its angular spin frequency. 

Combining advantag es of the outer gap and TPC models, 
lOiao et all (2004. 2007) originally suggested the annu l ar gap 
mode l , whic h has been further developed by IDu et alj d2010t 
1201 ll 120121) . The site for generation of high energy pho- 
tons is mainly located in the vicinity of the null charge sur- 
fac ffi Being cons istent with the physically calculated spec- 
tra (Du et al. 201 1), the Gaussian emissivities are numerically 
assumed when simulating light curves. The key emission ge- 
ometry parameters a and £ are not convincingly confirmed so 

4 They are defined as a set of special field lines that satisfy the condition 
of SI ■ B = at the light cylinder. 

5 It is defined as a surface w here the Goldreich- Julian (GJ) charge density 
Pgj 1 Goldreich & Julian 1 969 ) is zero. 
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FIG. 1 . — The calculated acceleration electric field on an open field line in 
the annular gap are plotted for typical MSPs with four groups of parameters 
of surface magnetic field B and spin period P. The electric field can extend 
from the pulsar surface to the null charge surface or even beyond it. The 
maximum electric field in the annular gap is sufficiently high to accelerate 
primary particles to ultra-relativistic energies. 

far, we adopt these two values from related literature if they 
exist. If there are not any reliable values, we just use some 
random values according to theories of MSPs' magnetic field 
evolution. By hypothesizing reasonable emissivities and the 
magnetic inclination angle, this model can make wide fan-like 
emission beams, thus well reproduce the relevant light curves 
cut by a suitable viewing angle. 

2.2. Acceleration Electric Field 

Co-rotating charge-separate d plasma is fi lled in 

a pulsar magnetosp h ere dGoldreich & Julianl 119691: 
IRuderman & Sutherland! 119751) . When reaching some 
regions near the light cylinder, charged particles can not 
exceed the speed of light (c), thus escape from the magneto- 
sphere. To compensate the escaping particles, the pulsar has 
to supply sufficient charged particles to its magnetosphere. 
This dynamic process continuously happens, thus a huge 
acceleration electric field is generated in the magnetosphere. 
This is the general mechanism for accelerati on electric field 
(£11 ), which is suitable for both young pulsars (IDu et al.ll201 ll 
120121) and MSPs. The charged particles with opposite signs 
are simultaneously exporting from the annular gap and core 
gap, and satisfies the condition of circuit closure in the whole 
magnetosphere. 

We assume the flowing-out particles' c harge density equals 
to the local Goldreich- Julian density dGoldreich & Julianl 
119691) at a radial distance of r ~ R^c, here R^c = j- is a ra- 
dius of light cylinder. T he detailed calcu lation method and 
formulae are presented in IDu et all (1201 ll) . The acceleration 
electric fields of typical MSPs are shown in Figure Q] using 
four sets of surface magnetic field B and spin period P. For 
MSPs with large or small values of P and B, we found that, 
in all cases, the electric field in the inner region of the annular 
gap is sufficiently high (E\\ > 10 6 Vcm" 1 ). 
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The charged particles accelerated in the annular gap or core 
gap are flowing out along a field line in a quasi-steady state. 
Using the derived acceleration electric field, we can obtain 
the Lorentz factor T p of primary particles from the balance of 
acceleration and curvature radiation reaction 

r p = (^)i=2.36xl0VX6 5 ' 

where e is the charge of an electron, p-j the curvature radius in 
units of 10 7 cm and Eng the acceleration electric field in units 

of K^VcnT 1 . 

The primary particles are accelerated to ultra-relativistic en- 
ergy with typical Lorentz factors of T p ~ 10 6 - 10 7 , because of 
the huge acceleration electric field in the annular gap. Since 
a lot of 7-ray photons are generated by the primary particles 
via curvature radiation and inverse Compton scattering pro- 
cesses, abundant e^ 1 pairs are subsequently created through 
two-photon annihilation and photon-magnetic-absorption (7- 
B) processes. 

3. SIMULATING RADIO AND 7-RAY LIGHT CURVES 
FOR MSPS 

Thanks to the ferm/-LAT, we now know that some MSPs 
are muti-wavelength emitters which have detectable radio and 
7-ray pulsed emission. According to the obs ervations of 
phase lag (A) betw een radio peak and 7-ray peak dAbdo et al.l 
1201 Oat iKerr et alj 120121) . MSPs can be divided into four 
classes. PSR J0034-0534 represents a class of MSPs which 
has phase-aligned light curves (A ~ 0); PSR J0101-6422 
stands for another class which has moderate radio phase lag 
A <~ 0.2-0.3 with quite complex radio or 7-ray light curves; 
PSR J0437-4715 stands for a third class which has larger ra- 
dio lag A ~ 0.43 and PSR J1744-1134 is a fou rth class of 
MSPs whose 7-ray peak precedes the radio peak ( Abdo et al. 
I20l0ah . we will model this MSP in future when high signal- 
to-noise 7-ray light curves are available. We process Fermi 
Pass 7 data to derive the observed light curves for three MSPs 
according to the radio timing solutions of MSPs from Fermi 
Science Support Center (FSSCfl We select events with ener- 
gies of > 0. 1 GeV within 2° of each MSP's position and with 
zenith angles smaller than 105°. The key filter conditions for 
good time interval are rock angle < 52° and angsep(RAMSP, 
DECmsp, RAsun, DEC sun ) < 5°, where RA and DEC are 
right as cension and declin ation respectively. Then we use 
tempo2 dHobbs et alj|2006h with Fermi plug-in to obtain the 
spin phase for each photon. Finally we obtain the high signal- 
to-noise 7-ray light curves for the three MSPs (see red lines 
in Figure 2, 3, 4). 

A convincing model should have simple clear emission ge- 
ometric picture with reasonable input parameters, which can 
not only reproduce multi-wavelength light curves for young 
pulsars but also for MSPs. In this paper, we will jointly 
simulate radio and 7-ray light curves for PSR J0034-0534, 
PSR J0101-6422 and PSR J0437-4715. We briefly introduce 
the simulation method here. As shown in table 1, a and £ 
are "first-rank" parameters, which are primarily adopted from 
the literature if they exist. When there are not any convinc- 
ing values, we tend to use reasonable values of a from rele- 
vant t heory on magnetic field evolution of pulsars (iRudermanl 
11991b and £ is adopted randomly according to the simulated 
emission pattern (photon sky-map). When a and £ are fixed, 

6 http://fermi.gsfc.nasa.gov/ssc/data/access/lat/ephems/ 



several other model parameters are carefully adjusted for the 
emission regions until the observed light curve of the corre- 
sponding band can be reproduced. The model parameters for 
three MSPs are listed in table 1 . 

3.1. PSR J0034-0534 

PSR J003 4-0534 is the ninth 7-ray MSP detected by the 
Fermi-LPK dAbdo et al.ll2010b1) . and it has strong 7-ray and 
radio pulsed emission with phase-aligned light curves. To re- 
veal the emission region of this MSP, we use the annular gap 
model to jointly model the radio and 7-ray light cu rves. The 
simula tion method is the same as described in §3. 1 o flDu et al] 
d201 11) . The key idea is to project all radiation intensities 
in both the annular gap and core gap to the "non-rotating" 
sky, considering some physical effects (e.g. aberration effect 
and retardation effect). Here we use numerical emissivities to 
speed up the calculations, since the assumed emissivities are 
consistent with the physica lly calculated spectra, as shown in 
Figure 8 of IDu et al.l d20Tlb . 

From Figure 2, we find that both radio and 7-ray emission 
are mainly generated in the annular gap region co-located at 
intermediate altitudes r ~ 0.24-0.56j?T c, whi ch leads to the 
phase-aligned light curves. lAbdo et alj d2010bl) used the TPC 
and outer gap geometric models with a = 30° and ( = 70° to 
obtain the light curves for PSR J0034-0534, and they derived 
similar conclusions that radio emission region extends from 
0.6/?lc to 0.8/?lc and 7-ray region extends from 0.12/?lc to 
0.9/?lc- It is found that this MSP has a lar ger transverse emis- 
sion region for radio emission. Moreover. lVenter et alj d2012l) 
developed the traditional outer gap and TPC model, adopt- 
ing the similar idea of numerically assumed emissivity of 
piecewise-function, and derived the phase-aligned radio and 
7-ray light curves for three MSPs including PSR J0034-0534. 
It seems that PSR J0034-05 34 has off-peak pulsed 7 -ray emis- 
sion up to 100% duty cycle dAckermann et alJ201 11) . which is 
not reproduced by the annular gap model. This might be due 
to the lack of knowledge on emission geometry and magnetic 
field configuration, and we will further develop our model to 
study this MSP in detail in the future. 

3.2. PSR J0101 -6422 

The observed light curves of PSR J0101-6422 have com- 
plex features: the 7-ray profile is likely to have three peaks; 
while the radio profile contains three peaks that occupy nearly 
the whole rotation phase (see left panels of Figure 3). To well 
model both light curves, we consider both cases of single-pole 
and two-pole emission-picture with many attempts on a large 
number of parameter space. The light curves of this MSP fa- 
vor a single-pole emission picture with a = 30° and £ = 48°, 
and the results are shown in Figure 3. For a viewing angle of 
C = 48°, the radio interpeak originates from the core gap with 
high altitudes of 0.3/?lc to 0.7^lc; while the other two peaks 
with a bridge (at the phases of ~ -0.3 and ~ 0.45) originate 
from the annular gap with altitudes of OASR^c to 0.78/?lc- 
The 7-ray profile is similar, the interpeak (at phase of ~ 0) 
mainly originates from the core gap with high altitudes of 
0.32/?lc to 0.6/?lc, while the other two peaks originate from 
the annular gap with altitudes of 0.08/?lc to 0.4/?lc- 

According to the annular gap model, we note that the ra- 
dio emission from PSR J0101-6422 is quite asymmetric in 
magnetic azimuth. This is possibly due to the special physi- 
cal coherence condition and propagation effects in the pulsar 
magnetosphere. We add some discussions on radio emission 
in §4. 
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FIG. 2. — Simulated radio and 7-ray light curves for PSR J0034-0534. The observations are presented in the left panels, and the radio data is taken from the 
Fermi ephemeris website http://fermi.gsfc.nasa.gov/ssc/data/access/lat/ephems/ The photon sky-map (middle panels) for an inclination angle of a = 35° and the 
corresponding simulated light curves (thin black hnes in the right panels) for a viewing angle of f = 70° are also presented, using the single-pole annular gap 
model that well reproduce the observed phase-aligned light curves (thick red lines in the right panels). 



TABLE 1 

Model parameters of multi-wavelength light curves for three MSPs 



Band 




A e 


fA 


o"e.A 




C0,C1 


CT e.c2 






J0034-0534 


(a = 


35°; C = 


70°) 






>0.1GeV 


0.75 


0.85 0.8 


0.32 


0.007 


0.15 


0.006 


0.0052 


Radio 




0.58 




0.013 




0.00065 


0.00053 






J0101-6422 


(a = 


30°; C = 


48°) 






>0.1GeV 


0.50 


0.85 0.8 


0.3 


0.009 


0.15 


0.006 


0.006 


Radio 




1.5 




0.006 




0.0075 


0.0078 






J0437-4715 


(a = 


25°; C = 


42°) 






>0.1GeV 


0.1 


0.1 1.3 


0.03 


0.002 


0.02 


0.001 


0.001 


Radio 




1.2 




0.0004 




0.0002 


0.0002 



NOTE. — a and f are magnetic inclination angles and viewing angles; k and A are two geometry parameters to determine the peak altitude in the annular 
gap; e is a parameter for the peak altitude in the core gap; <ta and <jc are length scales for the emission region on each open field line in the annular gap and the 
core gap in units of Rlc, respectively; cre.A is the transverse bunch scale for field lines in the annular gap; ag.ci and cr$ ci are th e bunch scale fo r field lines of 
-180° < ip s < 90° and 90° < tp s < 180° in the core gap, respectively. The detailed description of these symbols can be found inlDuet al. 2011). 



3.3. PSRJ0437-4715 

PSR J0437-4715 i s a very close MSP wit h a distance of 
0.16kpc to the Earth ( Manchester et al.ll2005l) . and has multi- 
wavelength emission. IChen et all ( I 19981) suggested that PSR 
J0437-4715 was an aligned rotator. We therefore adopted 
a relatively small magnetic inclination angle as a = 25° 
and a reasonable viewing angle as C = 42° from the high- 
precision radio timing observations (Ivan Straten et al.ll200lt 
iHotan et al.12006 ). Bogd anov et al.l (120071) also used the value 
of £ = 42° to successfully model the thermal X-ray pulsations 
of this MSP. We apply our annular gap model to jointly simu- 
late its radio and 7-ray light curves, and the results are shown 
in Figure 4. We emphasize that the radio interpulse can be re- 
produced by our model, although it does not precisely match 
the observations. The radio emission including the main peak 
and interpulse originate from a much higher and narrower re- 
gion in the annular gap region with high altitudes of 0.48/?lc 
to 0.57/?lc; while the 7-ray emission is generated in the an- 
nular gap region with lower altitudes of 0.0647?lc to 0.15/?lc 
located in the same magnetic pole. This leads to a large radio 
lag of A ~ 0.43. 



4. CONCLUSIONS AND DISCUSSIONS 

Pulsed 7-ray emission from MSPs has been observed by the 
sensitive Fermi-LKY. Particularly, the specific pattern of radio 
and 7-ray emission from the PSR J0101-6422 challenges the 
outer gap and TPC models. A convincing model should ap- 
ply not only to young pulsars but also to MSPs. In this paper, 
we used the annular gap model to jointly model the radio and 
7-ray light curves for three representative MSPs PSR J0034- 
0534, PSR J0101-6422 and PSR J0437-4715 with distinct ra- 
dio phase lags. For PSR J0034-0534 with phase-aligned radio 
and 7-ray light curves, both bands are mainly generated in the 
annular gap region co-located at intermediate altitudes. For 
PSR J0101-6422 with complex radio and 7-ray pulse profiles, 
we presented the best simulated results for this type of MSPs. 
The radio interpulse originate from the core gap at higher alti- 
tudes; while the other two radio peaks with a bridge originate 
from the annular gap region. The interpeak originates from 
the core gap region and the other two peak from the annu- 
lar gap region. For PSR J0437-4715 with a large radio lag, 
the radio emission (including the interpulse) originates from 
a much higher and narrower region in the annular gap region, 



Annular Gap Model for MSPs 



>0.1 GeV (Fermi) 



A A 





Pulsar Phase # 
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FIG. 3. — Similar as Figure 2, but for PSR J0101-6422. The radio data is obtained from Kerr et al. (2012). The inclination angle a = 30° and viewing angle 
£ = 48° are used to model the complex radio and 7-ray light curves with moderate radio lag for this distinct MSP. 
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FIG. 4. — Similar as Figure 2, but for PSR J0437-4715. From right bottom panel, the radio interpulse at phase ~ —0.4 and main peak at phase ~ 0.4 are 
reproduced by the annular gap model, and they are generated from a higher and narrower regi on in the annular gap region of the same magnetic pole as the 7-ray 
emission. The radio profile for this MSP is observed from Kunming 40-meter radio telescope (Hao et al. 2010). The inclination angle a = 25° and viewing angle 
£ = 42° are used. 



and the 7-ray emission has lower altitudes. 

From simulations of these MSPs, the annular gap model fa- 
vors a single-pole emission pattern with small inclination an- 
gles (a < 35°) for MSPs. This result is compatible with the- 
ories of magnetic field evolution of MSPs in binaries: some 
recycled pulsars tend to have aligning magn etic filed moment , 
i.e. small magn et ic inclinatio n angle a dRudermanl 1 1 99 It 
IChen et alj 11998b . lLamb et al.l d2009b presented a concrete 
discussion on the a evolution of MSPs while they were recy- 
cling in low mass X-ray binaries, and they note that the strong 
interactions between spinning superfiuid neutrons and magne- 
tized superconducting protons in a pulsar's core force the spin 
axis to change. A MSP would be an aligned rotator (a ~ 0°) 
if the star's north and south magnetic poles are forced toward 
opposite spin poles by the accretion disk, or would be an or- 



thogonal rotator (a ~ 90°) if both of the star's magnetic poles 
are forced toward the same spin pole. Moreover, it is certainly 
uncle ar what happens wh en a MSP's recycling process fin- 
ishes, lYoung et al.l ( 12010b analyzed the new pulse width data 
of normal pulsars, and found that the spin and magnetic axes 
would align when they spin down due to dipole radiation and 
particle outflowing. 

Several MSPs are simply assumed to be nearly orthog- 
onal rotators because they have a radio interpulse sepa- 
rated by a large p hase of >, 180 ° from its main pulse 
( IChen et all 1 1998b . iGuillemot etail d2012bb studied multi- 
wavelength light curves for two MSPs (PSR B 1937+21 and 
PSR B 1957+20), and found that fitting the radio polarization 
data of PSR B 1937+21 favors an orthogonal rotator. As with 
most RVM fits the confidence area is large but the orthogo- 
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nal configuration is further supported by the altitude-limited 
TPC and outer gap models. However, this is not universally 
true, at least in the case of the annular gap model. The radio 
light curve of MSPs (e.g. PSR J0437-4715) with an interpulse 
can be well explained by the annular gap model with a small 
magnetic inclination angle. 

By simulating light curves for MSPs in the annular gap 
model, we found that the radio emission mainly originates 
from the high-altitude narrow region in the annular gap re- 
gion. The radio emission pattern (photon sky-map) is patch- 
like in our model. The radio emissivities on each field line (in 
the annular gap or core gap) vary slightly (nearly uniform), 
but the case for 7-ray light-curve simulation is quite different, 
they vary a lot. High energy (7-ray and X-ray) emission is 
generated by non-coherent radiation from relativistic primary 
particles and pairs, while radio emission is suggested to be 
generated by coherent radiati on due to two-stream instabilit y 
of outward and inward pairs dRuderman & Sutherlandl[T975l) . 
Here we focus on studying radio and 7-ray light curves for 
MSPs, the concrete emission mechanism including polariza- 
tion, spectral properties and long-term stabilities of radio lags 
are however needed to further study, considering the c oher- 
ent condition and propagation effects. lHan et alTd 19981) sys- 
tematically studied the radio circular polarization for pulsar 
integrated pulse profiles, and found that sense reversals of cir- 
cular polarization are observed across the conal emission in 
some cases, unrestricted to core components. The polariza- 
tion property of high-altitude radio emission from both annu- 



lar and core gaps is a valuable subject to be investigated in 
future. We will keep on improving our model to present bet- 
ter simulated light curves, especially for the phases of leading 
wing of peak 1, trailing wing of peak 2 and off-peak pulses. 

In this paper, we simulated radio and gamma-ray light 
curves for 3 MSPs. When high signal-to-noise data at other 
wavelengths is available in the future, we will re-simulate 
light curves and fit the multi-wavelength phase-resolved spec- 
tra. In sum, the annular gap model is a self-consist ent model 
not only for young pulsars dDu et al.l 1201 ll 120121) . but also 
for MSPs, and multi-wavelength light curves can be well ex- 
plained by this model. 
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